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SUMMARY 


A rapid scanning two-dimensional laser velocimeter (LV) has been used to measure 
simultaneously the vortex vertical and axial velocity distributions in the Langley Vortex 
Research Facility. This system utilized a two-dimensional Bragg cell for removing 
flow-direction ambiguity by translating the optical frequency for each velocity component, 
which was separated by band-pass filters. A rotational scan mechanism provided an 
incremental rapid scan to compensate for the large displacement of the vortex with time. 
The data were processed with a digital counter and an on-line minicomputer. Vaporized 
kerosene (0.5-pm to 5-jim particle sizes) was used for flow visualization and LV scat- 
tering centers. The overall measured mean-velocity uncertainty is less than 2 percent. 
These measurements were obtained from ensemble averaging of individual realizations. 

INTRODUCTION 

The laser Doppler velocimeter (LV) has rapidly become an accepted diagnostic 
tool in the field of fluid dynamics research. There are basically two types (i.e., the ref- 
erence beam (ref. 1) and the dual scatter or "fringe" system (ref. 2)). The fringe system 
has been used more extensively because of its superior performance capability and rela- 
tive ease of application as compared with the reference beam system. It exhibits better 
signal-to-noise ratios and can be operated either in the backscatter mode or the forward 
scatter mode. (The geometry of the test environment often dictates the operating mode 
best suited for a particular measurement situation.) Some of the important features of 
the LV are as follows: 

(1) It is a nonintrusive technique which is less likely to disturb the flow than mea- 
urement devices such as hot-wire probes. 

(2) The velocity is measured directly (i.e., independent of other flow parameters 
such as temperature and pressure), which consequently eliminates the need for complex 
instrument calibrations or the need for other quantitative parameters. 

(3) Good spatial resolution can be obtained for detailed flow-structure studies. 

(4) The LV system is highly responsive to local-flow velocity fluctuations in both 
magnitude and direction. 

As the measurement requirements in complex aerodynamic flow have increased, so 
has the design sophistication of LV instrumentation to comply with these needs. For 



instance, during the past several years, vortex research has been expanding rapidly 
(refs. 3 to 11). Most of this research expansion has been attributed to the advent of the 
wide-body jet transport for which lift-induced vortices generated by the aircraft present 
a potential hazard to other encountering aircraft. This phenomenon has led to the need 
for instrumentation to measure the vortex velocity structures and to obtain quantitative 
information on the vortex characteristics. A number of methods have been investigated. 
For example, probes such as hot wires (refs. 12 to 14) and pressure sensors (ref. 15) 
have been used. These devices, as previously stated, are subject to producing flow dis- 
turbances and often require calibration before the data can be interpreted. This is par- 
ticularly true in more severe test environments where the meandering motion of the 
vortex system can further complicate the data reduction. 

The applications of LV to swirling flows have been documented. Lennert et al. 

(ref. 16) used a nonscanning LV system for mapping the flow field of wing-tip vortices of 
a flying aircraft. A scanning LV system, approximately 1 scan/sec and translation 
speed/scan up to 1.6 m/sec, has been used for measuring vortex velocities in a wind 
tunnel (ref. 17) and a water towing facility (ref. 18). However, in a static-air towing 
facility where the vortex drift (downwash) is relatively large and produces transient type 
measurements, much higher scan rates are needed. To accurately define the vortex 
velocity distribution in this type facility, scan rates on the order of 30 Hz (scans/sec) 
were required. This report describes a rapid scanning two-dimensional (i.e., vertical 
and axial) LV system developed for the Langley Vortex Research Facility and the related 
systems' tests and evaluations. Typical results of these tests are given. 

A discussion of the tracking fidelity of the LV scattering centers (particles) in the 
vortex flow is presented in the appendix written by Michael J. Walsh. 

Use of trade names or names of manufacturers in this report does not constitute an 
official endorsement of such products or manufacturers, either expressed or implied, by 
the National Aeronautics and Space Administration. 

SYMBOLS 


A aspect ratio 

b model wing span, m 

C D drag coefficient 

lift coefficient 
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D 


particle diameter, m 


D y sample volume diameter, m 

fg effective Bragg cell modulation frequency, MHz 

f c counter reference clock frequency, MHz 

f^ Doppler frequency, MHz 

V = *B ± f d 

M number of threshold crossings 

m exponent 

N total number of individual velocity measurements 

N Re r relative Reynolds number 

n refraction index 

r vortex radius, m 

S distance of apparent shift of focus, m 

fringe spacing, mm 
T r counter reset time, sec 

t window thickness, mm 

u half -cone angle of laser beam, deg 

Up angle of laser beam to the normal to the window, deg 

V particle velocity, m/sec 

Vj r fringe velocity, m/sec 
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V- individual measured velocity, m/sec 

model velocity, m/sec 
V mean velocity, m/sec 

V m measured velocity, m/sec 

x distance in direction of center line of vortex, m 

y spanwise distance (coordinate in crossflow plane of vortex), m 

y c lateral distance of vortex core from plane of symmetry, m 

z coordinate in crossflow plane of vortex 

a angle of attack, deg 

/3 angle of incidence, deg 

/ 3 1 angle of refraction, deg 

y turbulence intensity 

9 angle between intersecting laser beams, deg 

X laser wavelength, m 

p, viscosity, N-sec/m^ 

p density, kg/m^ 

a standard deviation 

Subscripts: 
g gas 

p particle 
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X 


x-component 


y y- component 

z z -component 

An arrow over a symbol represents a vector quantity. 


VELOCITY MEASUREMENT TECHNIQUE 

The velocity measurement technique may be better understood by considering the 
basic principles of the laser velocimeter, which depends on Mie scattering and the 
Doppler effect. Shown in figure 1 is a simplified illustration of the fringe model where 
two parallel coherent beams of radiation of equal intensity are focused by a lens to form 
"interference fringes” at the focal plane (sample volume). The fringe pattern is in the 
form of equal light and dark regions. The fringe spacing Sf is given by 

(1 

2 sin (0/2) v 


where X is the wavelength of the laser radiation and 0 is the angle between the inter- 
secting beams. 

Individual particles (typically several micrometers in diameter) embedded in the 
moving fluid pass through the fringes, thus causing light to be scattered. This scattered 
light is detected by a photodetector, thereby generating a fluctuating electrical current 
in the form of a burst signal at a periodic frequency f^ given by 



where |v| is the absolute normal velocity component of the particle crossing the fringe 
pattern. The Doppler shift principle can also be used to describe the LV measurement 
technique (ref. 19), where the frequency of the incident light from each of the beams 
scattered by the particles in the flow is changed by the movement of particles through the 
beams. By mixing these frequencies optically, a difference frequency proportional to the 
velocity of the particle can be detected, thus yielding the same results as shown in 
equation (2). 

Velocity directionality can be obtained with the use of a Bragg cell. This device 
produces moving fringes, with a velocity proportional to the Bragg cell modulation fre- 
quency. The Bragg Doppler-shifted frequency f^' is given by 
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(3) 


where fg is the Bragg cell modulation frequency and f d is the Doppler frequency. 
Particles moving in the same direction as the fringes will produce a frequency shift 
lower than the modulating frequency. The converse is true for particles moving in the 
opposite direction of the fringes. 


Optical System 

The optical system used in the present tests (fig. 2) is a confocal system configured 
to operate in the backscatter mode (i.e., utilizes a common lens arrangement for both 
transmitting and receiving). It is similar to that of reference 20, with the significant 
difference being a unique incremental scanning mechanism (ref. 21). 

The laser velocimeter uses an argon-ion laser that operates in the TEMqq mode 
and provides up to 5 W of continuous power at the 514.5-nm wavelength. The laser output 
beam is focused by a 2 -diopter lens into the center of a two-dimensional ultrasonic Bragg 
cell. The Bragg cell allows two orthogonal velocity components (axial and vertical) to be 
measured. The cell diffracts the incoming laser beam into an array of several diffraction 
orders. The diffracted orders of the laser beam have the Bragg cell frequencies added 
to the laser's normal optical frequency. The zero order, (0,0), has only the optical fre- 
quency, whereas the diffracted first orders, (0,1) and (1,0), have the optical frequency 
plus the Bragg cell frequency for that axis. The focusing of the laser beam by the 
2-diopter lens minimizes the laser-beam diameter near the origin of the angular diver- 
gence of the diffracted orders. This results in a minimum beam diameter for the pro- 
jected focus in the test section. The Bragg cell is positioned so that the maximum energy 
is diffracted into the shifted orders. The power to the drivers is then adjusted until the 
optical power in the (0,0) order, the (0,1) order, and the (1,0) order is equal. The power 
in the cross order, (1,1), is then also equal. The equalized intensity in the laser beams 
produces maximum fringe contrast at the projected focus (sample volume). A more 
elaborate discussion of this procedure is given in reference 22. The Bragg cell is driven 
at 15 and 25 MHz for the axial and vertical components, respectively. These frequencies 
were selected so that each velocity component could be retrieved from the photodetector 
signal by band-pass filters with minimum cross talk. The four orders from the Bragg 
cell are then collimated by a 1.5-diopter lens, transmitted through a four-hole detector 
mirror, and imaged into the dead air space of an optical cell by a 20-diopter FI: 1.5 lens. 
The optical cell is used to eliminate extraneous Doppler signals from particles passing 
through this primary focus. The beams emerging from the optical cell then pass through 
the scan wheel. This wheel is a 40.5-cm-diameter disk containing 16 windows. The 
windows vary in thickness from 0 to 30 mm, in 2 -mm increments. All windows are 
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fabricated from schlieren-grade optical crown from the same melt, with both faces 
parallel to 5 arc sec of wedge or less, flat to 1/10 wave of mercury green light, and 
antireflection-coated to reduce the effects of aberrations. 

The astigmatic and spherical aberrations introduced by the windows can be calcu- 
lated by knowing the cone angle of the beam, the thickness, and index of refraction of the 
window (ref. 23). (See sketch (a).) 



The spherical aberration (S.A.) is given by 


t_ 

n 



n cos u 




sin^ u , 


(4) 


where u equals half of the cone angle of the beam. For the system discussed herein 
with a cone angle of approximately 1°, an index of refraction of 1.57, and a window thick- 
ness of 30 mm, the spherical aberration is approximately 1.7 pm. The astigmatism is 
given by 

-t 

y'n 2 - sin 2 u p 

where Up is the angle of the beam to the normal to the window which is, again, of 
thickness t and index of refraction n. For this system with the angle of the beam to 
the normal to the window of approximately 3.5°, the astigmatism is approximately 
41.8 pm. 


9 9 

n^ cos^ Up 

- sin ^ u. 


-1 


(5) 
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The scan wheel provides scan rates (longitudinal displacement of the focus along 
the optical axis) up to 30 times/sec, with 16 discrete individual measurement positions 
for each scan. It operates on the refraction principle, whereby the insertion of plane 
parallel windows (of various thicknesses) between the primary focus in the optical cell 
and a projection lens arrangement produces an apparent incremental longitudinal dis- 
placement of the primary focus. The displacement due to refraction is shown in 
sketch (b). 



Sketch (b) 


If the thickness of the window is t, the refractive index n, the angle of incidence ( 3 , 
and the angle of refraction / 3 ' , then the displacement S may be shown to be 


y = t tan ( 3 ' - (t - S) tan (3 


( 6 ) 


S = 


tan ( 3 ' 
tan /3 


and by Snell's law 


S 



cos <3 
cos /3’ 


( 7 ) 


( 8 ) 


This displacement is magnified (by a factor of approximately 50) by the 2 -lens projection 
system. The first projection lens is a 100-mm-diameter FI: 4 lens corrected to better 
than \/8 wave-front deformation. The second projection lens is an FI: 5.6, with an 
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effective focal length of 127 cm. The two lenses are separated by about 76.2 cm. This 
separation can be varied to allow the positioning of the projected LV focus across the 
test section. The scan range encompassed by the higher speed scan is 0.4 m at a nomi- 
nal 3-m projected distance. The sample volume formed in the test section is typically 
0.2 mm in diameter and 1.4 cm long at the 3-dB level. The actual sensing-sample 
volume depends on the discrimination level of the processing instrumentation. The 
sample -volume length given was measured by moving a diffuse scatterer through the 
sample volume and normalizing the detected output. The sample-volume diameter was 
measured with a microscope. Particles in the test section crossing the projected focus 
(sampling volume) scatter light which is collected and imaged back through the lens 
system and scan wheel to emerge as a collimated beam from the FI: 1.5 lens. The col- 
limated beam is directed by the flat detector mirror to a 190-mm focal-length lens. 

This lens images the collimated beam onto a,105-/im-diameter pinhole. The light 
passing through the pinhole is detected by a photomultiplier. 

Signal Conditioning Electronics 

The signal conditioning electronics shown in figure 3 consist of a photodetection 
assembly, a digital frequency counter, and a minicomputer plus peripheral equipment. 

A high-gain, low-noise photomultiplier with a S-20 wavelength response was used 
to detect the light scattered by the particles passing through the sample volume and col- 
lected by the optical system. The photomultiplier output was coupled to a matching 
50-ohm line-driver amplifier. 

Band-pass filters (25 MHz and 15 MHz) were used to separate the vertical and 
axial velocity information transmitted by the Bragg cell carrier frequencies. Mixers 
were employed to downshift the frequency spectrum for each component by driving them 
at 20 MHz to produce an effective Bragg frequency of 5 MHz. This was done to reduce 
the probable operating error of the digital frequency counter. 

A simplified diagram of the digital counter-type processor is shown in figure 4. 

It is a time -domain device which uses a direct counting approach for the time it takes a 
particle to traverse a preset number of successive fringes. This approach eliminates 
most of the data ambiguity derived from frequency domain processors such as spectrum 
broadening effects because of the finite transit time for a particle to pass through the 
sample volume (ref. 1). Other advantages include the lower input signal-to-noise ratio 
and duty cycle (Doppler burst rate) required, as opposed to spectrum analyzers and fre- 
quency trackers. The input signal, in a form of a Doppler burst, is filtered to remove 
the pedestal voltage. The resulting signal is applied to a double-threshold zero-crossing 
detector, a dual-function device used as a wave shaper and noise rejector for reducing the 
likelihood of false data. A simplified operation of this network is as follows: The first 
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zero crossing is detected only after the input signal with a positive slope crosses the 
negative threshold level. In a similar fashion, the second zero crossing is detected 
when the input signal with a negative slope exceeds the positive threshold level. Con- 
tinuing this process generates a series of square pulses of constant amplitude. The 
first two pulses reset and arm the clock-counting circuits. The third pulse starts the 
period-counting process. Additional noise rejection is achieved with a data validation 
circuit. It has two identical timing circuits (i.e., a 5- and 8-cycle time count) driven 
with a 500-MHz reference clock (2-nsec resolution). The output for each circuit is 
applied to a comparator network. If the resulting 5:8 comparison falls within the pre- 
selected limits (0 to 10 percent), the data are accepted as valid. 

The validated data go to a storage register which has a binary digital output and is 
interfaced to a minicomputer. There are 10 data lines and 4 exponent lines for covering 
a Doppler frequency range of 1.2 kHz to 100 MHz. Additional lines for external inhibit 
and synchronization are also available. The frequency f^’ is determined from the 
conversion equation 

„ 3.2 x 10 10 

Frequency = 

2 m x reading 

where the values of the exponent m and the reading are indicated by the digital counter 
processor. These values go into the minicomputer (binary form), which is programmed 
to calculate the frequency from the conversion equation. 

A 4-bit optical encoder mounted on the scan wheel is used to identify each measure- 
ment location. The output from this device, along with the output data from the digital 
counters, is fed into a 16-bit minicomputer which catalogs the related velocity measure- 
ments in the flow. The data-acquisition operation of the minicomputer is initiated by a 
photoelectric trip switch activated by the passage of the strut support for the vortex- 
generating model. The trip switch is located with the plane of the measurement region. 
After the switch is activated, the computer remains inhibited until all the output lines 
from the 4-bit encoder are low (indicating the first measurement location). This time 
delay is stored in the computer. The total time is calculated from scan speed and delay 
time and is used to correlate measured velocity with vortex age. The mean velocities 
and associated standard deviation can be computed, tabulated, and plotted as a function of 
time or position. 


SYSTEM PRECISION 

The instrumentation error introduced by the laser velocimeter depends basically on 
the fringe spacing f(X,0) and the measured frequency. The laser is equipped with an 
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etalon in the cavity space to provide a stable excitation-line width of less than 0.001 cm -1 . 
Therefore the effects of fringe spacing due to the uncertainty of the wavelength X is 
extremely small and is therefore neglected. The cross-beam angle 0 for each mea- 
surement location was determined within 0.007°. The velocity bias error due to cross- 
beam angle is expressed as 


^M = |0,5cot (0/2) A0 


( 9 ) 


A typical cross-beam angle range from 0 = 1° to 0.5° and A 0 = 1.15 x 10 4 rad yields 
a velocity error of approximately 1 percent. 

The velocity error associated with the frequency determination arose from the 
signal-processing electronics. The error contributed by the signal-processing elec- 
tronics (digital counters) depends on the quality of the input signal and the internal control 
parameters of the counter (refs. 24 and 25). For the measurements presented herein, the 
signal level was generally better than 15 dB above the noise. This level was sufficient to 
set the counter threshold for discriminating against the noise background. 

The main internal control parameters of the counter that affect the system preci- 
sion are digital quantization and clock synchronization. The following assumptions and 
procedures outlined in reference 25 are used to approximate the bias and random error 
contributed by this device. Digital quantizing error (random) occurs during the measure- 
ment of the particle transit time (which consists of discrete clock pulses) through a given 
number of fringes. The synchronization error (bias) develops during the start of the 
timing sequence where the reference clock frequency and the signal frequency are 
matched. These error sources in the measured Doppler frequency f^ can be 
expressed as follows: 

Quantization error — 


where 


V m (f) - V(f) 

W) 



( 10 ) 


5 = 


1+i 


( 11 ) 


and 
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( 12 ) 


d_ 

Mf, 


Clock synchronization error - 


v m (f) - v(f) 
v(f) 



(13) 


Terms in equations (10) to (13) are: f^, true mean Doppler frequency; fg, effective 
Bragg cell frequency; f c , counter reference-clock frequency; M, number of input signal 
cycles used for processing; y , flow turbulence intensity; V m , measured velocity; and 
V, true velocity. Based on the previous equations, as well as an average Doppler fre- 
quency of 0.5 MHz, and effective Bragg cell frequency of 5 MHz, 20-percent flow turbu- 
lence intensity, 8 -input signal cycles, and 500-MHz reference clock frequency, estimates 
of these errors are made and summarized in the following table: 


Source 

Bias, 

percent 

Random, 

percent 

Cross-beam angle 

±1.0 


Quantization 


Negligible («0.01) 

Clock synchronization 

.0687 



APPLICATION TO VORTEX FLOW FIELD 

In order to evaluate the system performance in a swirling flow environment, tests 
were conducted in the Langley Vortex Research Facility. A schematic diagram of this 
facility is shown in figure 5. This facility uses a gasoline-powered carriage to tow air- 
craft models through a static -air test area at a constant velocity. The overhead track 
for the carriage is 549 m long. The test area is covered (except a 5 -cm opening for the 
model support) to provide isolation of the test-model-generated vortices from the wake 
of the carriage. The test section is 91.37 m long, 5.49 m wide, and 5.18 m high. For a 
more detailed description of this facility, see reference 8. The test model (fig. 6) used 
for generating the wake vortices is a 0.03-scale model of the Boeing 747. It has a wing 
span of 1.79 m and is equipped with leading- and trailing-edge flaps to simulate a landing 
or cruise configuration. 

The seeding material used in these tests was vaporized kerosene, which has a mass 
density of 0.82 g/fn^ and a size distribution of 0.5 pm to 5 pm, as shown in figure 7. The 
size distribution was determined with a commercial par tide -sizing system. 
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The following test procedure was used to obtain a velocity profile measurement of 
the wake vortex. The vortex-generating model is towed through the test section at a con- 
stant nominal velocity of approximately 30 m/sec. (See fig. 8.) As the generated vor- 
tex system drifts through the scan region at a nearly constant rate, it is rapidly scanned 
(30 times/sec) over a distance of approximately 0.4 m. The rapid scan is needed because 
of the time dependency of the vortex characteristics. Measurements were made in 
regions covering approximately 1 semispan. Although velocity measurements were 
made for the entire run time, only vortex core penetrations are presented. Typical mea- 
surements at any location during a run were based on approximately 20 velocity deter- 
minations during a 2.0-msec observation time. The separation distance of the model 
from the vortex measuring region was varied by changing the height of the model and 
keeping the height of the sample regions constant (2.3 m). The model heights above the 
sample regions were preselected for the purpose of obtaining measurements in the near 
field (roll-up region) and the far field (along the stable region). These tests were con- 
ducted with the model at representative lift coefficients corresponding to the model test 
configuration. For instance, tests for the landing (flaps) and the cruise (basic aircraft) 
configurations were with lift coefficients of 1.2 and 0.67, respectively. At the test loca- 
tion, heated kerosene vapor was deployed for flow-field visualization (fig. 9) and LV scat- 
tering centers. The resulting flow-visualization data as a function of time are recorded 
on videotape. A video timer (also initiated by the photoelectric switch) was used for the 
time display. These measurements complemented the LV measurement by pinpointing 
the time that the vortex core penetrated the sample region, whereby the corresponding 
scan number could be determined. Using this technique greatly reduced the data search 
time. 


VELOCITY MEASUREMENT ACCURACY 

The operating principle of the LV technique depends partly on Mie scattering and 
involves measuring the velocity of particles suspended in the moving fluid in order to 
obtain the velocity of the fluid. In determining the overall measurement accuracy based 
on this concept, consideration must be given to the particle-flow tracking fidelity. An 
analysis of particle dynamics in a swirling or vortex flow is discussed in the appendix. 

It was found that particles (<1 pm) in the vortex core migrate in the radial direction, 
which results in the lack of data for far -field velocity measurement in the core region. 
The resulting conclusions of this analysis imply that less than 1 -percent error is intro- 
duced in the final measurement due to particle -tracking error. Therefore the overall 
uncertainty in the measured mean velocity, including particle tracking, cross-beam angle 
measurements, and frequency determination errors, is less than 2 percent. 
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The measurement data from the digital counters represent the detection of the 
occurrence of individual events in the flow, whereby the final velocity measurement 
result is based on an arithmetic ensemble averages. A given set of these individual 
realizations is expressed as 


N 



i=l 


a = 



(15) 


with the statistical uncertainty in the mean given as ± ~^=. Data acquired by means of 

v/n 

individual realizations, which are subject to biasing effects (i.e., randomly occurring 
velocity-dependent events in an assumed uniformly seeded flow) will cause the resulting 
velocity measurement to be biased toward the higher velocities. On the other hand, in 
flow situations where Bragg cells are used, the velocity measurement is biased toward 
the lower velocities. Thus, the combined effect has a tendency to reduce the overall 
bias. There are many proposed treatments for the velocity biasing problem (refs. 26 
to 30). McLaughlin and Tiederman (ref. 26) have shown that the corrected velocity 
(unbiased) can be obtained from weighted ensemble averages, as shown in equation (16) 



N 

i 

N 

I 

i-1 


i w i v i 


1 w i 


(16) 


where the weighting function W i is expressed as 1 / |Vj . Meyers' equation (ref. 27), 
which takes into account the Bragg frequency, is given by 


V = 


N 

1 w i v i 

i=l 

N 

1 w i 
i=l 


(17) 
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where 


w 10 S f + T r (v fr + Vj) 
i_ D v (v £r + Vi) 

Sf fringe spacing 

T r counter reset time 

Vf measured velocity 

Vf r velocity of fringes 


(18) 


D y sample volume diameter 

Corrected velocity data were computed based on equation (17), which takes into account 
the effects of the Bragg cell. These calculated values were compared with the arithmetic 
ensemble averages, and the results revealed insignificant differences. Therefore the 
data presented in this report are derived from arithmetic ensemble averages. 


RESULTS AND DISCUSSION 


Wake-vortex flow velocity profile measurements were made at a model speed of 
30 m/sec. Plots for the cruise and landing configurations are shown in figures 10 and 11, 
respectively. The measured vortex velocity distribution for the cruise condition was 
obtained for a 0.67 lift coefficient and separation distance (between the model and the 
measurement plane) of 7.1 spans. In a similar fashion, measurements were made for 
the landing configuration, where the lift coefficient is 1.2 and the separation distance is 
55 spans. The lack of data points in the vortex core is due to the absence of sufficient 
scattering centers in this region. Particles are centrifuged out of the core as the flow 
develops, with the larger particles (>2 (im) migrating out first. For the cruise condition, 
the vortex drift (downwash) was low enough that several successive scans for the same 
run could be used for mapping the velocity distribution for essentially the same spatial 
location. However, for the landing condition where the vortex displacement with time 
was much larger, separate runs were required. The theoretical curve, superimposed on 
the measured values, is based on the equation, given in reference 31 as 


Voo 7T 2 A 


y -y c 

( 6b ^ 

1/2 + b 

|y - y c |< 

0 y - y cl j 

y - y c 


(19) 
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where C L is the lift coefficient, A is the aspect ratio, b is the model span, y is 
the spanwise distance, and y c is the lateral distance of the vortex core from the plane 
of symmetry. This equation is the same as used in reference 31 to analyze data obtained 
with a similar model in a towing basin, and the measured velocity data were in good 
agreement with the theoretical values over the region compared. 

During the simulated landing configuration, several vortices were shed from each 
wing of the test model. The near-field (x/b ~ 1) vertical velocity distribution of the 
dominant vortex, produced by the outboard flaps, is shown in figure 12. Typical vortex 
vertical velocity distributions are shown in figure 13 at various separation distances 
between the model and the measurement region for the basic aircraft. These measure- 
ments were made at a nominal lift coefficient of 0.67. The simultaneously measured 
vertical and horizontal velocity distributions of the wake-vortex system are shown in 
figure 14. These values are normalized to the model velocity. The peak velocities for 
the vertical and horizontal velocity components are approximately 25 and 10 percent of 

the model velocity, respectively. The statistical uncertainty of the mean -£=■ for both 

\JN 

velocity distributions, when normalized to the model velocity, occurred within the plotting 
symbols. 


CONCLUSIONS 

A rapid scanning two-dimensional laser velocimeter (LV) system has been devel- 
oped and tested in the Langley Vortex Research Facility. On the basis of the results 
obtained, the following significant observations can be made: 

1. The major instrumentation error source was the determination of the cross- 
beam angles, d- v which contributed approximately 1-percent bias error to the system; 
radiation wavelength and electronic counter contributions were negligible. 

2. The overall measurement accuracy involving particle dynamics, measurement 
statistics, and velocity bias for particles with size distributions of 1 to 3 pm and a 
specific gravity of 0.82 g/cm^ contributed less than 1 -percent error due to flow tracking 
fidelity; and the data treatment using ensemble averages with the application of Bragg 
cells did not require data correction due to velocity biasing effect. 

3. The LV system is capable of simultaneously measuring both vertical and axial 
flow velocity components in a near field or far field of the vortex system, whereby quan- 
titative information can be retrieved. 

4. Particles (>1 pm) in the vortex core migrate in the radial direction which results 
in the lack of data for far -field velocity measurement in the core region. 
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5. The rapid scan rate minimized the effects of the vortex relative high-drift rate 
(downwash) so that an accurate definition of the vortex system could be obtained. 

6. The measured velocity data were in good agreement with the theoretical calcu- 
lated values over the region compared. 

Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
March 19, 1980 
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APPENDIX 


PARTICLE MOTION IN A VORTEX 

Michael J. Walsh 
NASA Langley Research Center 


The following section analyzes the three-dimensional motion of various size par- 
ticles passing through a vortex in air. 

The general governing equation given by Soo (ref. 32) for a particle traveling in air 
reduces to the following equation for the special case of a spherical particle traveling in 
air: 


dV 


c^n t 


^ = 3 -D-Re ? r ^g /- 
dt 4 p D 2 l g P 
P P 


where the relative Reynolds number N Re r is defined as 


p V - V D 

N = S- S Pl_J2 

Re ’ r 


(Al) 


(A2) 


Equation (Al) assumes that viscous drag is the only force acting on the particle. As dis- 
cussed in reference 33, equation (Al) can be treated in component form as 


dV n.x 3 C D%e.r M V g .x ~ V o.x) 
dt 4 p d 2 


P P 


(A3) 


dV. 


P»y _ 3 C D N Re,r ^g( V g,y ~ V p,y) 


dt 


p D *■ 
P P 


(A4) 



dt 


3 C D N Re,r ^g( V g,z 


P„D. 


P P 



(A5) 


For the present analysis, x is in the direction of the vortex axis, and y and z are in 
the cross-flow plane of the vortex. To complete the system of six differential equations 
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that are needed to solve for the three-dimensional particle motion, three additional equa- 
tions are obtained from the definition of velocity 


dx = y 
dt p,x 


(A6) 


& = y 
dt p,y 


(A7) 


dz = y 
dt P,z 


(A8) 


Thus, the three-dimensional particle motion through a vortex can be computed by a 
numerical procedure using the gas flow-field velocity V„, viscosity p density p , 
particle density p p and size D^, the drag coefficient C D , and the initial conditions of 
particle position and velocity. For the present calculations, the gas velocity distribution 
in the cross-flow plane of the vortex is the theoretical one shown in figure 11 as a function 
of radial position from the center of the vortex. The velocity component parallel to the 
vortex axis is equal to 30.9 m/sec. The gas density and viscosity were evaluated at 
atmospheric pressure and a temperature of 300 K. The particle density was equal to 
919.5 kg/m3. The drag coefficient was evaluated from the following equation given in 
reference 34: 



24 

Re,r 


^1 +0.1 5N 


0.687^ 
Re,r j 


(A9) 


which gives good predictions of incompressible sphere drag for r < 200. For the 

present calculations, the particles were initially at zero velocity prior to being accel- 
erated by the vortex at x/b = 0. The calculations indicated that velocity lags were 
small (generally less than 1 percent for 1- to 3-p.m particles, even though the particles 
migrate in the radial direction. This means that accurate laser velocimeter (LV) mea- 
surements can be obtained with 1- to 3-/im particles. 

As previously mentioned, the 1- to 3-pm particles do migrate in the radial direc- 
tion. It is of interest to see if the particle-motion calculations explain the lack of 
LV data for r/b < 0.02 at x/b = 55, which is the measurement station for the LV data 
presented in figure 11. Changes in the particle's radial distance from the center of the 
vortex are presented in figure 15 for various size particles passing through the vortex. 
The data indicate that the change in radial position depends on the particle size and initial 
radius at x/b = 0.0. The dashed line in figure 15 represents the minimum r/b at 
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which LV measurements could generally be obtained during the present tests. Note that 
the data in figure 15 are limited to x/b < 6 because the calculations were too time- 
consuming to continue them all the way to x/b = 55. It was believed that enough infor- 
mation could be obtained by extrapolating the r/b values for x/b < 6 to x/b = 55. 

The calculations indicated that all the particles greater than 2 pm migrate out of the 
region r/b< 0.02 by x/b = 55. The only particles left in the region r/b < 0.02 at 
x/b = 55 are 1 pm or smaller than were originally at r/b < 0.01. This represents a 
reduction of a factor of 4 in the number density of 1-pm particles for r/b < 0.02. Thus, 
the particle -motion calculations confirm the results obtained in the LV measurements. 
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Figure 1.- Laser velocimeter technique. 
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Figure 2.- Optical system. 
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Figure 3.- Data acquisition and processing system. (PMT is photomultiplier tube 

and BP is band pass). 
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(c) Digital pulse train from zero-crossing detector. 


Figure 4.- Digital counter-type processor. 
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Figure 5.- Langley Vortex Research Facility. 





Figure 6.- The 0.03-scale model. 
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Figure 7. 
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Figure 9.- Flow visualization of vortex system. 





Normalized vertical velocity 



Figure 10.- Vertical velocity distribution at x/b = 7.1, 
a = 6.33°, Cj^ = 0.65, and = 30 m/sec (cruise 
configuration). 
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Normalized vertical velocity. 


Fairing 



y/b 

Figure 11.- Vertical velocity distribution at x/b = 55, a = 4°, 
C L = 1.2, and = 30.9 m/sec (landing configuration). 
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Figure 14.- Vertical and axial velocity 
center at x/b = 30; o; = 3.5°; C L 
configuration). 


39 






1. Report No. 2. Government Accession No. 

NASA TP-1661 

3. Recipient's Catalog No. 

4. Title and Subtitle 

A SCANNING LASER-VELOCIMETER TECHNIQUE FOR 
MEASURING TWO-DIMENSIONAL WAKE-VORTEX 
VELOCITY DISTRIBUTIONS 

5. Report Date 

May 1980 

6. Performing Organization Code 

7. Author(s) 

Luther R. Gartrell and David B. Rhodes 

8. Performing Organization Report No. 

L-13426 

10. Work Unit No. 

505-31-23-05 

9. Performing Organization Name and Address 

NASA Langley Research Center 
Hampton, VA 23665 

11. Contract or Grant No. 

13. Type of Report and Period Covered 
Technical Paper 

12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, DC 20546 

14, Sponsoring Agency Code 

15. Supplementary Notes 


Appendix by. Michael J. Walsh. 

16. Abstract 

A rapid scanning two-dimensional laser velocimeter (LV) has been used to measure simul- 
taneously the vortex vertical and axial velocity distributions in the Langley Vortex Research 
Facility. This system utilized a two-dimensional Bragg cell for removing flow-direction 
ambiguity by translating the optical frequency for each velocity component, which was sepa- 
rated by band-pass filters. A rotational scan mechanism provided an incremental rapid 
scan to compensate for the large displacement of the vortex with time. The data were 
processed with a digital counter and an on-line minicomputer. Vaporized kerosene 
(0. 5-fim to 5-pm particle sizes) was used for flow visualization and LV scattering centers. 
The overall measured mean-velocity uncertainty is less than 2 percent. These measure- 
ments were obtained from ensemble averaging of individual realizations. 


18. Distribution Statement 

Unclassified - Unlimited 

Subject Category 35 

page) 21. No. of Pages 22. Price* 

40 $4.50 


17. Key Words (Suggested by Author (s)) 

Laser velocimeter 
Velocity measurements 
Test techniques 
Wakes 
Vortices 

20. Security Classif. (of this 

Unclassified 


19. Security Classif. (of this report) 
Unclassified 


For sale by the National Technical Information Service, Springfield, Virginia 221S1 


NASA-Langley, 1980 


















National Aeronautics and 
Space Administration 


THIRD-CLASS BULK RATE 


Washington, D.C. 
20546 

Official Business 

Penalty for Private Use, $300 


Postage and Fees Paid 
National Aeronautics and 
Space Administration 
NASA-451 





POSTMASTER: 


If Undeliverable (Section 158 
Postal Manual) Do Not Return 




